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We consider ideal fluid and equivalent scalar field dark energy universes where all four known 
types of finite-time, future singularities occur at some parameter values. It is demonstrated that 
pressure/energy density of such quintessence/phantom dark energy diverges in spherically-symmetric 
spacetime at finite radius or at the center. This may cause the instability of the relativistic star or 
black hole in such universe. The resolution of the problem via the extra modification of the equation 
- - - of state is briefiy discussed. 
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^ \ I. INTRODUCTION 

■ The discovery of the late-time universe acceleration brought to the playground the number of dark energy (DE) 
models with the effective equation of state (EoS) parameter w being very close to —1 in accordance with observational 
data. It is known that phantom/quintessence models lead to the violation of all/some of the energy conditions. Such 
models unHke to ACDM with w — ~1 lead to the number of quite surprisingconsequences in the remote future. For 
instance, phantom DEs are characterized by the future Big Rip singularity Some of quintessence DEs bring the 
O ^ universe to softer finite-time singularity in the future. Such finite-time future singularities may represent so-called 
1^ ' sudden singularities 0,01 or some other singularity types which are classified in ref.[J|. It is evident that the presence 
^ . of finite-time future sinffularity in the course of the universe evolution may show up at the current epoch. One example 
' has been given in ref.pf where it was conjectured that sudden singularity^ of specific modified gravity DE may make 

, the relativistic star formation process being unstable. The resolution of the problem is to introduce the higher-order 
J> ' curvature terms [3| relevant only at the early universe in such a way that future singularity disappears. 
^-H ■ In the present work we consider the specific dark fiuid which contains all four known finite-time singularity types 
. The reformulation of it as scalar DE model with the same FRW asymptotic solutions for the corresponding scalar 
' potentials is also made. The energy density /pressure of such singular DE may become divergent in the spherically- 
; . symmetric spacetime at finite radius or at the center. In a sense, that is the way the finite-time singularity manifests 
CO ' itself as singularity of spherically-symmetric space. This may lead to the instability of relativistic stars (in the same 
, way as for modified gravity DE model in ref.[§]) or instability of black holes located in such dark energy universe. 
It indicates that number of current DEs with such properties may be problematic for realistic description of current 
accelerating universe. Some extra EoS modification by the terms relevant at the very early universe may be necessary 
in order to resolve this problem. Such modification is discussed briefiy in the last section. The reconstruction method 
to find the specific dark energy responsible for any singularity of spherically-symmetric space is also presented. 



I 



0\ 

o 



X 



II. FINITE-TIME SINGULARITIES IN THE DARK FLUID UNIVERSE WITH THE EXPLICIT 

EQUATION OF STATE 

In the present section we show the appearance of the finite-time singularities in the illustrative but sufficiently 
realistic dark fiuid with the following equation of state (EoS) : 

P^~P + Ap^, (1) 

with constant A and a. We work in the spatially fiat FRW space-time 

ds^ ^ -dt^ + a{tf {dx'f . (2) 

i=l,2,3 
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Then since H = a/a, the conservation law in the FRW universe, p + 3H {p + p) = 0, gives 3^4^ = which can be 
integrated as 



(3) 



Here Oq is a constant of integration. In case a = 1, the EoS |(T|) reduces to the usual perfect fluid with constant EoS 
parameter w = —1 + A. The universe is regular in the future when w > —1 but there appears the Big Rip singularity 
when w < —1 for the constant EoS parameter dark fluid. 
The solution of the FRW equation 



is given by 



when a = 1 , A > 
when a = 1 , ^ < 



ao (to - t) ^^^^ 

When a ^ 1, by defining a new variable N, which is called e-folding, the FRW equation ([4]) can be rewritten as 
Then in case 1/2(1 — a) = 1 or a = 1/2, we find 



(4) 



(5) 



(6) 
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On the other hand, when a ^ 1/2, 



Then the scale factor a is given by 



a = age 



2(1 -g) ^^^4^ _ VSKAjt-tp) 
I -2a 2 



aoe«— , C^jh-^[^y-{\-a)^{-iA)^ 

3(1-°) _ / ,\Ti^ 2(l-e») 

oraoe^(*o-*) 1-- ' / ,n , 



, C. 
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For the first case in ([9]) , one may put io = without loss of the generality. We may summarize the behavior of a as 
follows. 



oo (^0 - ty'^^^^ , 



a = < 



ooe 



Nqc 1 — 

2(1-°) 
aCt 1-2° 



C: 



2(1-°) 
C(to-t) 1-2° 



or aoe 

Then the Hubble rate is given by 



^(^)-^(i-.)^(-3A)- 



when a = 1 , ^ > 
when a = 1 , A <Q 

when a = i , A <Q 
when a ^ 1, i 



C^-T^ f (a-i)— (-3^)-- 



(10) 



H = I 
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2(l-a) ^,l/(l_2a) 
1-2q 



when a = 1 , ^ > 
when a = 1 , ^ < 
when a = ^ , A <Q 

when a ^ 1, i 
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Now one can describe the future, finite-time singularities of the universe filled with above dark fiuid. When a < 0, 
there occurs Type II or sudden future singularity [3, When < a < 1/2 and 1/(1 — 2a) is not an integer, there 
occurs Type IV singularity. When a = 0, there is no any singularity. When l/2<a<lora = l and A < 0, there 
appears Type I or Big Rip type singularity. When a > 1, there occurs Type III singularity. 
The above general classification of singularities was first given in 

• Type I ("Big Rip") : For t ^ tg, a ^ co, p ^ oo and \p\ — > oo. This also includes the case of p, p being finite at 

• Type II ("sudden") : For t ^ ts, a ^ as, p ^ ps and \p\ oo 

• Type III : For t ^ ts, a ^ as, p co and \p\ — > oo 

• Type IV : For t — > t^, a — > Os, p 0, |p| and higher derivatives of H diverge. This also includes the case 
in which p (p) or both of p and p tend to some finite values, while higher derivatives of H diverge. 



In case of Type II singularity, where a < 0, H vanishes as ~ {to — t)^^^'^ when t ^ to and therefore we find 



Here, ts, as{^ 0) and ps are constants. 

In case of Type II singularity, where 
p vanishes as it follows from the FRW equation Then, near the singularity, the EoS |(T]) is reduced to 

p ^ Ap" . (12) 



On the other hand, in case of Type I singularity, where l/2<a<lorQ; = l and A < 0, H and therefore p diverge 
when t ^ to- Then the EoS ^ reduces to 

p -p or p -(1 - yl)p . (13) 

In case of Type III singularity, where a > 1, H and p diverge when t to and therefore the EoS |[T]) reduces to 

P^Ap<- . (14) 

Hence, we illustrated the occurrence of all four types of future singularity for the above dark fiuid. Actuallj^, the 
corresponding singular asymptotics of such as well as of the more complicated EoS dark fiuid are given in ref. 



III. FINITE-TIME SINGULARITIES IN THE SCALAR FIELD DARK ENERGY UNIVERSE 

We now consider the scalar-tensor theory, which also leads to the above singularities. One starts from the following 
action for the scalar-tensor theory; 

s = jd\\^-^^-'^d,4>d^4>-V{<t>)] . (15) 

Here u;(0) and V{(j)) are some functions of the scalar field 4>. Note that scalar field may be always redefined so 
that kinetic function is absorbed. In fact, in case w((/)) is positive definite, if we redefine the scalar field <j) by 
= / d(l)\/\uj{4>)\, we obtain 

^ = / rf'^ { ^ T \d,vd^^ - ^(^) } . (16) 

Here the potential V{if) is defined by V{lp) = V {(f) {f)). The minus (plus) sign in =F corresponds to positive (negative) 
Let us consider the theory in which V^(0) and are given by 

^(0) = , V{cty) = \ [3/(^)' + /'(</.)] , (17) 

where /((/>) is a proper function of (f). Then the following solution is found [H] 

= t , H{t) = fit) . (18) 
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For the action (flSl) . the Hubble rate H in (fTTj) could be realized by the following scalar models: 



_3A 



a;(0) 



4(l-a)C</.]^ 
l-2a 



1/(0) 



3Aj 
4 



27^^^ _ 3A 



Noe- 



when a — I 



when a = - , /I < 
2 



6(1 - a)C^(j)l-'- + C^l-'° 



when a ^ 1, 



(19) 
(20) 

(21) 



Here </>+ = and <j>- ^ to - 0. In ((191), ('^-) when ^ > (>1 < 0). 

Let us investigate the form of V{ip) in l|16p when we rewrite the action lfT5|l . 



27A^ 
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3C^ 



when a = 1 



when a = i 
when a = I 



when a 



(22) 
(23) 
(24) 

(25) 



Here 4>+ — 4> and (f>- ^ to — if), again, 4>o is the constant of the integration, the upper (lower) sign corresponds to — 
(+) sign in ([16]). In ([23]), and the upper (lower) sign corresponds to A > (A < 0), iVo < (TVq > 0), 
C > (C < 0), C > ((7 < 0), respectively. 
We also find when a ^ 1,^, |, 



4^iF(l-2a)(l-a)C 



3 -4a 
(1 - 2af ^2 



3-4q 

, 2-4a 



4WiF(l-2a)(l-a)C 
^ = ^ 7Z (*o 



6 (1 - a) 



(3 - 4a) 



4^T(l-2a)(l-a)C 



(3 - 4a) (p 



4^T(l-2a)(l-a)C7 



,(26) 



3-4o 
1 2-4a 



T/(^) 




(1 - 2a)^ K 
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(3 - 4a) if 



C 



(3 - 4a) (y5 



4A/=F(l-2a)(l-a)C, 



4^iF (1 - 2a) (1 - a) C 



4(l-a) 



(27) 



C > 



The upper (lower) sign corresponds to — (+) sign in (fT6|) and — i_2a ^ > — 1-2° ^ '^j (|26|l . 
(-^^?<0) in {22!). 

Note that the such scalar theory exactly reproduces the singularity and/or behavior in the scale factor lfTO|) and the 
Hubble rate ifTTj). which is generated by the EoS H]) of previous section dark fluid. This shows that such a perfect fluid 
with the EOS ^ can be realized by the scalar field with specific potential. The important lesson of this presentation 
is that big class of ideal fluid as well as scalar fleld dark energies brings the future universe to the flnite-time singularity 
of one of the four types under consideration. 



IV. SINGULARITIES OF SPHERICALLY-SYMMETRIC SPACETIME FILLED WITH DARK ENERGY 

In this section we show that finite-time singularities of dark energy models in above two sections manifest themselves 
as radius singularities of spherically-symmetric spacetime filled with such dark energies. In it has been pointed 
that curvature singularity is realized inside the relativistic star for a viable class of /(i?) -gravities (for review of 
viable, realistic models of that sort, see Q)- Since such spherically-symmetric solution with a naked singularity is 
inconsistent, this result indicates that large star (or even planet) could not be formed, or such a relativistic star could 
be unstable in such a theory. 
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It is known that viable modified gravity also shows all four above types of finite-time singularity Hence, it is 
natural to expect that qualitatively similar situation should be typical for any dark energy which brings the universe 
to finite-time singularity. Motivated with these observations, we investigate the EoS dark fiuid, which generates 
curvature singularity in spherically-symmetric spacetime. Especially in this section, we investigate the singularity, 
which is generated for a finite (non- vanishing) value of the radius. The singularity at the origin (vanishing radius) is 
investigated in the next section. 

Let us first consider what kind of (perfect) fiuid could generate a singularity. We concentrate on the singularity 
which occurs for a finite radius for spherically symmetric solution. Assume the metric has the following form: 

ds^ = -e-^-^de -f e^'-^^Ur^ + r^d^l . (28) 

Here d^\ expresses the metric of two-dimensional sphere. Eq. ((28|) expresses the arbitrary spherically symmetric and 
static space-time also in the presence of some matter, that corresponds to the inside of the star or planet. Then the 
Einstein equations have the following form: 

1 j\ „2A 1 

K^e^^ , (29) 
^Ve^^ , (30) 

/• ar I— 
d^i/ _^ f d'l^ dX\ /du 2 2X 

dr"^ \dr dr J \dr r J 

Here p is the energy density and Pr and Pa are the radial and angular components of the pressure. In the following, 
as for the usual perfect fiuid, we assume Pr = Pa- If one does not impose this assumption, we can consider general 
types of singularity. By combining ((30l) and ([31]) and deleting p = p,- = Pa, it follows 

d^v f dv\^ dv f \ dX\ 1 dX 



IdX 




r dr 




1 di^ 




r dr 










. dr r 



dr^ \dr ) ^ dr \r ^ dr J ^ r dr r^ ^ ^ 

The following kind of singularities at r = rg may be now considered: 

oo oo 

A(r) = Ao + Ai In h A„ (r - ro)"~^ , v{r) = vq + vi In + ^ Vn{r - tq)"""^ . (33) 

" n=2 " n=2 

By substituting (pS]) into (|32| . one obtains 

vi-vl + ViXi ^^^^ - 2iyiiy2 + 1^1X2 + 



[r - rof r-ro 



2Ai 



-1^1 -X2,y2 + - '- + - + O{{r-ro)) = 5 ) . (34 

Since the l.h.s. in l(34|) contains only the power (r — ro)™, where m is an integer greater than or equal to —2, 2Ai 
should be also an integer greater than or equal to —2. Furthermore if we assume 2Ai = —2, from the coefficients of 
(r — ro)~^, we find — h'f + e^^", which is inconsistent since the r.h.s. is positive definite. Therefore 2Ai must be an 
integer greater than or equal to —1; 2Ai > — 1. Then from the coefficients of (r — ro)~^, again, we find 

= 1^1 (1 — 1^1 + Ai) , that is, i^i = or t^i = Ai + 1 . (35) 

If 1^1 = Ai = 0, there is no singularity and we do not consider this case. In case I'l = and 2Ai = —1, from the 



coefficients of (r — tq) in l(34|) . we find ^ ^ ^ = ^- case i^i = and 2Ai > 1, one finds 

1^2 = -— ■ (36) 

The case ^ and i^i^Ai + 1 = 1/2 corresponds to the black hole, where tq corresponds to the horizon radius and 
it follows — 1!^2 + ^ ^ ^7^7" = 0. In case vi = \ and Ai = 0, we find ^ — 2v2 -I- A2 = 0. In case i^i = Ai + 1 > 3/2, 
one gets "^^^^ - 2vxV2 + 1^1X2 + 1^2X1 = 0. 
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We now investigate how p and p behave. In case vi = Q and 2Ai = —1, from Eqs. ((29l) and (|30|) . p and p are not 
singular and behave as p ~ — p ~ i^^h^- Therefore there could not be the singularity at r = tq. In case i^i = and 

2Ai = n > 1, Eqs. m and m give p ^ , P ^ ro(r-ro)" ' ^^^rc relation (136|) is used. Then p and p 

satisfy the following asymptotic equation of state (EoS) : 

P ~ -Cp^ . (37) 

Here C is a positive constant. Since p and p diverge at r = rg, there is a curvature singularity at r = rg. As we will 
see soon, the EoS ([ST)) corresponds to that in lfT4|) . which generates type III singularity. 

In case Ai = and vi — 1 case, we find p ^ i— + t- — ^— j— ) , P ^ i — — ^. Here p is finite although 

p diverges at r = tq. The divergence of p generates the curvature singularity. The EoS has the following form: 
pip — Po) ^ const.. 

In case 2Ai = n > and i^i = Ai + 1, we find 

K'^ (r — roj K [r — ro) 



and therefore the following asymptotic EoS: 



P={l + -]p. (39) 



n 



The divergence of p and p at r = tq means the curvature singularity. Since the EoS parameter w = p/ p = 1 + 2/n is 
positive, the EoS does not generate any finite-time singularity and corresponds to a = 1 and —1 + A—1 + 2/n. 

Hence, rather general case with the singularity occurence for finite r is investigated. We found here essentially two 
types of singularity expressed by the EoS fiuid l(37|) or l(39|) . These EoS fiuids could be compared with the asymptotic 
EoS dark fiuid generating the finite time singularity lfT2|) . (fT3| . and lfT4|) . First, one notices that EoS fiuid (|39| for 
finite radius singularity does not exist for dark fiuid generating the finite time singularity as in lfT2l) . (flSl) . and (fT4l) . 
(However, it could be that other dark fiuid generating future singularity shows up radius singularity in this example). 
Eq. l|12p has a similar structure but in (fT2l) . w = p/ p < —1 although w > 1 in l(39|) . The EoS (|37l) has a similar 
structure with those in lfT2|) and (fT2l) if we identify 

n 

a = . (40) 



In case of lfT2|) . a is negative but in case of (fT4l) . a is positive. Then the fiuid with EoS (fT4|) . which generates type 
III singularity, may generate finite radius singularity. Note, however, we do not have the explicit proof that the fiuid 
generating finite-time singularity always generate finite radius singularity, and vice versa. That depends also from the 
distance between energy scale which is typical for future singularity and the star/black hole formation energy scale. 
The important lesson of our study is that dark fiuid generating finite-time singularity may manifest it also as radius 
singularity of spherically-symmetric solution (i.e. stars, black holes, etc). 

In the above analysis, it is assumed that A(r) and i^(r) are real and therefore e^'*'^'') and e^''^''^ are positive. As we 
have only investigated the local behavior, the global structure of the spacetime is not so clear. Since e^^^''^ and e^"^^^ 
are positive, however, the number of the horizons must be even and in principle, the singularity can be naked. If there 
is a usual matter, whose distribution is also spherically symmetric, as a main component of the star, the singularity 
might be hidden inside the matter. If the radius of the matter becomes smaller than the radius of the singularity, 
however, due to the strong gravitational force near the singularity, the matter can shrink very rapidly, which may 
produce the strong radiation, say. X-ray. 

If e^'*'^''^ and/or e^''*^''^ are negative near the singularity, however, there must be horizon, where e^^^*"^ and/or e^'^'^''-* 
vanish, since e^^^''^ and e^''^''^ are positive in a region far from the star. We should note that if e^^° is negative, 
= 1/^+ e^^" has non-trivial solution even if 2Ai — —2. In case of l(38|) . where 2Ai = n > and i^i = Ai + 1, if e^^° is 
negative, p becomes negative and inconsistent. Therefore the singularity with 2Ai = n > and i^i = Ai -I- 1 could be 
forbidden in some cases (compare with singularities classification in ref.[^). 



V. CENTRAL SINGULARITY OF SPHERICALLY-SYMMETRIC SPACETIME FILLED BY DARK 

ENERGY 

In the previous section, we have investigated the singularity, which occurs at the finite (non-vanishing) radius. In 
this section, we consider the singularity at the center of the spherically-symmetric solution and clarify what kind of 
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EoS dark fluid could generate such a singularity. It is also shown that there are some exact solutions. Since such 
a solution with a naked singularity is inconsistent, even if the singularity exists at the center of the star, this result 
indicates that large star/planet could not be formed, or such a relativistic star is unstable. 

Assume that the energy density p and the pressure of the perfect fluid only depend on the radial coordinate r, then 
the conservation law gives 

dp dv , ^ , , 

0=^ + -b + p)^0, (41) 

which could be derived by using (|29| . l(30|) . and l|31|) . Especially for the perfect fluid with a constant EoS parameter 
w: p = wp, we flnd p — poc When is a monotonically increasing function of the radial coordinate r, as 

in the Schwarzschild metric, p becomes a monotonically decreasing function of r if w > or w < but p becomes 
a monotonically increasing function if > w > —1. Then the gravity effectively acts as a repulsive force when 
0>w>-l. 

For later convenience, it is convenient to write the EoS in the following form: 

p + p = F{p). (42) 
Here F{p) is an appropriate function. Then (|4T1) gives 

We should note that l(30|) with Pr = P = Pa can be rewritten as 

1 \ -n 

(44) 




Then by using l(29|) and (|44l) . one can delete A and obtain 



2 d^v 2_djy 2_ ^2dp 2_ (Idu , J_\ ]_\ -, 

^2p+l ^2p+l 



r dr 

By using l|42p and (|43l) . we may delete v and p in l(45 



-7f{p)dfi + 7FtpW~d:r{'^) +TW{pj^-^ 2 dp -2K=^g + 

~ _ 2 dp , 1 ^2„J.± F(r,)dr J. ± ' ^ ^ 



2 



d^p I 2 dF{p) fdp\ I 2 dp 2 



rF{p) dr 



Eq. l|46p can be regarded as the differential equation for p = p{r) for a given F(ji). With known solution of l(46|) with 
respect to p = p(r), the explicit form of v may be found by using l(43l) . Then by using (|44l) . we obtain A. Finally 
Eq. l(3T1) might determine some of the integration constants. 

As an example, we consider the perfect fluid with a constant equation of state (EoS) parameter w. It follows 
F{p) — ap, a= 1 + Then one of the exact solutions in l|46p could be given by 



which gives 



2r 2w r ^ ^ 1 / + Aa - 4\ 1 f + 6w + 1 

— in — In — , A = Ao = - m ^ = - in ^ — 

a ro w + 1 ro " 2 \ J 2 \ (w + if 



When r ~ 0, the scalar curvature R behaves as i? = — c^^'^" + ^) Hence, except a = —2, that is, w = —1/3, 
there is surely (naked) singularity at r = 0. 
We now consider more general case, as in ([33)) . 



A(r) =Ao + Ailn^ + ^A„r"-i , K^) = t^o + ^^i In - + i^„r"-^ (49) 



'^o :r2 ^0 
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By substituting (|49|) into (|32|) . one gets 



i/i - i/f + 1/1 + z/iAi + Ai + 1 , c2^"+^;^'=2^"(''^'-n)" ' /^rV^ 



+ 0(r-i)= -2 - ■ (50) 

Since the l.h.s. is only given by the integer power, in order that there could be a solution, 2Ai should be an integer 
n: 2Ai = n. We also find, if n < 0, there is no solution in l(50|) . therefore n > 0. 

When n = 0, by requiring the term with to vanish, we obtain e^"^" = (\/2 + 1 — i^i) [i^i + -\/2 — l). Since the 
r.h.s. is positive, -\/2 + 1 < i^i < \/2 + 1. When n = case, ([29]) and ([30]) gives 

(z.2_ 21^1 + 1) r2 ' K2_2t.i + l)r2 ' ^^^^ 



which gives the following EoS parameter. 



Unless ui = Q, p and/or p diverge at r = and therefore there is a curvature singularity there. As we require p 
could be non-negative, we find a constraint vi < Q ox vi > 2. Combining this constraint with the previous one 
-V2 + 1 < t^i < \/2 + 1, we find 

-\/2 + l<i/i<0or2<t/i<\/2 + l. (53) 

As w in l(52|) is constant, this case corresponds to a = 1 case in |(T]). Combining l[52|) with l(53|) . one gets | < if < 
2 + 3-\/2, — 3\/2 + 2 > w > — oo. In the former case, w is positive and therefore, there does not occur any finite-time 
singularity. In the latter case, however, w is negative and — 3-\/2 -I- 2 < — 1 and therefore this case corresponds to the 
Big Rip or Type I singularity. 

On the other hand, when n > 1, we find 



4 - n ± + 32 , , 

VI = (54) 



Then from ISSl) and dSOl), we obtain 



2 



Ai - 1 _„ / r \ " n — 2 ^ -» \ " '^i ~ 1 f ^ \ " %A? + 32 



which gives the following EoS parameter 



(55) 



1^1 -1 -n ± + 32 
"^a7~T= 2(n-2) ' 

which is constant again and corresponds to a = 1 case in Here Eq. l|54p is used. In order that p is positive, one 
finds n > 2. When n = 2, p vanishes and 



2c- 

p = , K^p ^ , ^ . (57) 

Since p is finite although p diverges, this case corresponds to Type II singularity, which corresponds to a < in the 
EoS H]). In the upper sign (-1— sign) in (|56| . w is positive and therefore this case does not correspond to any finite-time 
singularity. In the lower sign ( — sign) in ((56l) . however, w is negative and it corresponds to Type I singularity. 

Since p and/or p diverge at r = in above cases, there is a curvature singularity there. If usual matter is a main 
component of the star, the singularity could be hidden even if the singularity could be naked when there is no matter. 
When the usual matter is stripped, if the interaction between the usual matter and the anomalous matter, which 
could generate the singularity, is weak, the naked singularity could appear. Due to the strong gravitational force, the 
singularity could attract the matter around the singularity and strong radiation cold be generated. 

We have investigated rather general singularities. In order to treat further general cases, one can now consider a 
reconstruction program. When we require p behaves as 

p = P{t) , (58) 
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one can try to find F{p) which realizes (|58|) . Defining a function G{p) by 



by using ((58|) . Eg . (1461) can be rewritten as 

^ dGcPG K^r^P'{r) + 2 f dGV 
dr dr"^ r {K'^r'^P{r) + 1) \ dr J 

Eq. l|60p is the first order differential equation for dG/dr. By solving i(60|) . one may find G as a function of r: G — Gir). 
Then by using l(58|) . we get G as a function of p by G{p) = G (P~^(p)) . Here P^^{p) is an inverse function of p = P{r)- 
Using l(42l) and (|59l) . we find that the EoS is given by 

1 = (P + P) J" dp'G{p-\p')) . (61) 

Using the above equation, one may find the EOS dark fiuid generating general singularity expressed by fSS]) . 



VI. DISCUSSION 

In summary, we investigated dark fiuid/scalar DE universe where all four known types of future singularities occur 
for the corresponding choice of EoS parameters/scalar potentials. The energy density /pressure of such DE is singular 
in the spherical spacetime at finite radius or at the center. As a result, the relativistic stars and/or black holes in 
such universe may become unstable which makes the corresponding DE model to be problematic as realistic theory 
for current universe acceleration. One way to overcome this problem is to modify the initial EoS in such a way that 
finite-time singularity does not occur or occurs in far remote future so that it does not infiuence star/BH formation 
processes. 

The contributions to the FRW equations from usual matter could be neglected near the singularity. Then combining 
the FRW equations with l(42|) . we find —2H/k^ ^ F{p). Hence, if the absolute value oiF{p) is bounded, there does not 
occur Type I, II, or III singularities where H diverges. One can achieve this by extra modification of the corresponding 
EoS. From another side, in the same way as in second work of ref.[3| one can add inhomogeneous (gravitational) terms 
which are relevant at the early universe, to the EoS so that the future universe becomes regular. However, the 
introduction of such terms effectively means the modification of gravity (for general review of various modified gravity 
models see [0|). Eventually, this brings us back to modified gravity as consistent DE model. 
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